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^5 , We derive the dynamic and tiiermodynamic stability conditions for dilute multicomponent Bose-Einstein 

' condensates (BECs). These stability conditions, generalized for n-component BECs, are found to be equivalent 

' and are shown to be consistent with the phase diagrams of two- and three-component condensates that are 

, derived from energetic arguments. 

, ^ ■ PACS numbers: 

^ ■ I. INTRODUCTION 

m ; 

Interpenetrating superfluids have intrigued physicists for decades. The possibility of a superfluid mixture of two different 
isotopes of Helium motivated many early theoretical studies, including a hydrodynamic approach |1] and an approach using 
the Bogoliubov analysis \2^. (See also |3i] for a generalization of the Bogoliubov approach to mixtures in a different context.) 
"q , However, these mixtures have never been experimentally realized in superfluid Helium. 

^ ; Recent advances in creating and manipulating Bose-Einstein condensates in ultracold alkali gases have made it possible to 

d , probe the properties of interpenetrating superfluids by studying, for instance, two internal states of a single atom |4], different 
^ ' atomic species 011^, different isotopes of a single atom, or a combination thereof. One such advancement is the achievement 
I of a condensate of Ytterbium ff] which, with its five naturally occuiTing bosonic isotopes, opens the door for the creation of 
' O interpenetrating condensates with up to five components. Another is the development of optical Feshbach resonance techniques 
^ ] isl]. Unlike magnetic Feshbach resonance techniques, where typically only one scattering length can be tuned at a time, these 
new optical methods would, in principle, allow one to tune both the interspecies and intraspecies scattering lengths to probe a 
much more significant portion of parameter space |9]. 

In this paper, we investigate the fundamental stability and phase diagrams of such multicomponent dilute condensates, ex- 
eluding effects such as external potentials and metastability. We begin by deriving a simple relation to determine the thermo- 
dynamic and the dynamic stability conditions of an n-component miscible mixture, and we show that the stability conditions, 
I at least within our approximations, are equivalent. Then, through an energetic approach, we discuss these phase diagrams of 
two-component mixtures and those of three-component mixtures and demonstrate that these are consistent with the stability 
ff^ . conditions. 

o 

^ . 

O . II. STABILITY CONDITIONS 

a 

Let us consider a general miscible n-component condensate mixture in a box of volume V with periodic boundary conditions 
lIToll in the thermodynamic limit described by the following Hamiltonian 

O . i? = E / d^r -^t(j.)__^^(r) + ^^l{r)i^l(v)Mr)Mr) + E ffyV^J(r)^](r)^.(r)^,(r) (1) 

• H where, for simplicity, we assume that all species of atoms have the same mass m. We also assume that the number of particles 
rN of each species is fixed and cannot be converted into a different species. The atomic interactions are characterized by the 
^ interparticle contact pseudopotential where the coupling constants, gij = Anrh^ aij / m, can be given in terms of the two-body 
■ ■ ■ ' s-wave scattering length . 

Let us expand the bosonic field operators in terms of plane waves, ijji (r) = aj^.e*'^ '"/ W, where annihilates a particle 
of species i in state k and obeys the usual bosonic commutation relations. Let us assume that each species is condensed such that 
a macroscopic number of atoms are in state k = 0. Following the standard Bogoliubov approximation 1 1 1 ], we assume Oq and 
dg^ to be commuting c-numbers equal to \/N^, where N'l^ is the number of particles in the condensate of species i. Assuming 
that — Nq iV*, where = Nq + J^k^o ^-k '^k number of particles in species i (alternatively one could 

introduce a chemical potential in the grand canonical ensemble with given by miminimizing the thermodynamic potential). 
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we can approximate the Hamiltonian to leading order as 



(2) 



where pi = Nq/V. Note that the kinetic energy plays a negligible role because we consider only atoms in the k = state. The 
energy density of an n-component miscible mixture in the mean field approximation can therefore be written as 



Cmi ^ n ^ n 



(3) 



In this homogenous system, thermodynamic stability implies that the Hessian matrix defined by dp^ dp- is positive semi- 
definite m 



I.e. 



all eigenvalues of the matrix are non-negative, which is still possible in certain parameter regimes even if the 
interspecies scattering lengths are negative. Because dp-dp -E^^ = gij, thermodynamic stability of this rt-component miscible 
mixture implies that the matrix gij must be positive semi-definite. 

Now let us consider the dynamic stability of the system. Specifically, we investigate the conditions under which a small 
density perturbation increases exponentially in time or, equivalently, the excitation spectrum becomes imaginary. We derive the 
excitation dispersion relation by expanding the Hamiltonian of the system to the next order in the Bogoliubov approximation 
described above to arrive at 



H « Hq + i?intra + ^^intcr 

where the intraspecies contribution to the Bogoliubov Hamiltonian is given by 



N 



) + Pz9tz{a^a\ + aka!_k) 



= 1 k^O 



— {hkY I2m. The interspecies contribution is given by 

-Winter 



^ ^ V^Wgy 



V 



(4) 



(5) 



(6) 



In this approximation, the Hamiltonian is quadratic and, as such, can be diagonalized by canonical transformation. Therefore, 
one can transform this Hamiltonian describing weakly interacting particles into a Hamiltonian describing non-interacting quasi- 
particles with n independent modes of oscillation. The excitation energies hujn{k) can be derived by assuming a harmonic time 
dependence g^'^nC^)* of the particle operators. Writing the Heisenberg equations of motion for the particle annihilation and 
creation operators, we obtain 



hujn(k)a^ = (ifc + .gripOok + guPiali, + 2^ y ' « 



(7) 



The excitation energies are thus given by the eigenvalues of the 2n by 2n matrix given by the above equations and are given by 

[Sc^„(fc)]2 = iO.(<o+2A„) (9) 

where hujn{k) is the excitation energy and A„ are the eigenvalues of the matrix Ay = pjgij- If any of these eigenvalues 
are negative then the excitation energy becomes imaginary and the system becomes dynamically unstable in the sense that 
infinitesimal perturbations will grow exponentially with time. The dynamic instability starts at low momenta, and thus the 
condition for the system to remain dynamically stable is that the matrix Ay be positive semi-definite. Since pj are positive and 
positive semi-definiteness implies that all upper left submatrices are positive, the condition of Aij being positive semi-definite is 
equivalent to the condition of gij being positive semi-definite. 

The fact that the dynamic stability condition, which indicates whether a state is metastable or not, and the thermodynamic 
stability condition are the same implies that the energy density function's local and global extrema are the same. The energy 
surface is therefore downward convex for stable systems and upward convex for unstable systems 1 13]. The stability condition of 
this system, namely that g.^j must be positive semi-definite, has another important consequence. If the condition is met then the 
determinant of all upper left submatrices are positive implying that, for an n-component miscible system to be stable, all subsets 
of the species comprising the system must also independently satisfy the relevant stability conditions. A restrictive hierarchy of 
stability conditions is thus defined. 
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III. PHASE DIAGRAMS 



To illustrate these stability conditions, we will now discuss the ultimate fate of these mixtures when they are unstable; we 
consider one-, two- and three-component systems in that order We shall rely on purely energetic arguments in the mean-field 
approximation. We shall make the assumption that the system has been allowed a sufficiently long time to relax to its lowest 
energy state and ignore metastable states such as molecules. We ignore all effects from the thickness of the interface between 
regions and surface tension 1 17] as we assume that the kinetic energy plays a negligible role. We also assume that the coupling 
constants remain unchanged in the presence of other atoms whose effects we briefly discuss below. 

For the one-component case, the stability conditions yield the well known result that gn > must be positive in order for the 
condensate to be stable. If gn < then one can formally prove that the system will collapse in density space in a finite time if 
H{t = 0) < (where H is the full classical Hamiltonian) and in the thermodynamics limit 1 14]. In realistic systems, however, 
other effects such as three-body recombination will become important fTsil . 

For a system of two components, in this simplified mean field approximation we can safely make the assumption that only 
two stable states exist; the miscible case given by 

2 2 
cmi I 922P2 , 

^2 = — ^ + + .912P1P2, (10) 

and the completely immiscible case given by 

cim gilPl , 922pI r- ,,,, 

^2 = — ^ \ ^ h V5ll522PlP2- (11) 

This last equation can be easily derived if we assume that the total energy of the immiscible case is given by 

ET . ^ + ^ (12) 
' 2Va 2Vb 

where Va and Vb are the volumes of regions where species 1 and species 2 exist, respectively. We can simplify i?™ by noting 
that the total volume is fixed, i.e. Va + Vb — V, and, by either assuming that the system is mechanically stable Iil61 . i.e. 
dEy '^ / OVa = dE™ / OVb, or, equivalently, minimizing the energy with respect to Va under the constraint of Va + Vb = V 
IITIi . we arrive at eq. M \\ . 

It is clear that when > 5ii322 and gi2 is positive the system would energetically prefer to be in the immiscible phase. This 
condition is consistent with both stability conditions of the miscible phase derived above which indicate that if gjj > .91 1^22 
then the miscible system is unstable. If > 511522 and gi2 is negative then we would expect the system to collapse in a finite 
time because it can be shown that if i/(< = 0) < for a mixture, the system will exhibit finite time collapse ] 18]. 

We now turn to the rich phase diagram of systems with three interacting species. For simplicity we will assume that all species 
have an equal number of particles, i.e. Ni = N2 ^ N^, leading to a uniform constant particle density p. We also assume that the 
intraspecies coupling constants are normalized by the interspecies coupling constants and are equal, i.e. gn = 522 = 533 = 1- 
Here, the possible stable states are the miscible phase where energy density is given by 

Sr/P' = 3 + 2(512 +.923 +513). ^^3^ 
the completely immiscible phase where each species is separated from the other such that 

4'7p' - 1; (14) 

a phase in which species i is separated from a region where species j and k are miscible, with an energy density given by 

cim(») ; 2 _ 3 + 25j-fc + 2V2y/l + gjk 

^3 IP — 2 ' ^ ' 

and a phase where species i and j are miscible, i and k are miscible, but j and k are immiscible, with an energy density per atom 
given by 

. , . , , 3 + 2o„- + 2g,fc + + 
g^^Mip2 ^ £^ (jg) 

As in the earlier two-species example, we have assumed the pressures in each separate region to be equal, that is, we have 
assumed mechanical equilibrium. In the stable state, f™'-'''''', we have also assumed the chemical potentials of species i in the 



different regions to be the same. The phase diagram for the three-component case is shown in Fig. ^where we have assumed 
5i2 = 523 = 5 for simplicity. In this case, 5™'-^'^-' £ 
will collapse in density space if H{t = 0) < fisil . 



5i2 = 523 = 5 for simplicity. In this case, f™*^^'^-* and f ™'^''^'' are stable only when 513 = g. The shaded region in the figure 
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FIG. 1: The phase diagram for a three-component system is plotted against interspecies coupling constants gi2 = 323 ~ g and gis under 
the assumptions of Ni = N2 = N3 and gn = g22 = 533 = 1- Note that the shaded region is expected to collapse in density space if 
H{t = 0) < 0. 



Throughout this mean field analysis, we have assumed that coupling constants are determined entirely by two-body physics 
and remain unchanged in the presence of condensates. If many-body effects are taken into account, i.e. that the coupling 
constants may be modified in the presence of a foreign species, many exotic phases could occur, such as the possibility of a 
three-component condensate having a stable state with three distinct regions, each with different pairs of miscible species. In 
addition, other interesting effects would be possible, such as the presence of a third species making two otherwise immiscible 
species miscible, or the destabilization of a previously stable two species miscible mixture by the introduction of a third species. 
Finally, effects not considered in this paper, such as metastable states 1 19] and a non-uniform trapping potential |20], enrich the 
picture even further and come into play in realistic experiments involving trapped dilute condensates. 
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